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Abstract 



Results from the first study of isolated-photon+jet correlations in relativistic heavy 
ion collisions are reported. The analysis uses data from PbPb collisions at a centre-of- 
mass energy of 2.76 TeV per nucleon pair corresponding to an integrated luminosity 
of 150 jib^ 1 recorded by the CMS experiment at the LHC For events containing an 
isolated photon with transverse momentum pj > 60 GeV/c and an associated jet with 
Pj' > 30 GeV/c, the photon+jet pj imbalance is studied as a function of collision cen- 
trality and compared to pp data and PYTHIA calculations at the same collision energy. 
Using the p Z of the isolated photon as an estimate of the momentum of the associated 
parton at production, this measurement allows a characterisation of the in-medium 
parton energy loss. For more central PbPb collisions, a significant decrease in the 
ratio PjV Pj relative to that in the PYTHIA reference is observed. Furthermore, signif- 
icantly more p j > 60 GeV/ c photons in PbPb are observed not to have an associated 
p\ 1 > 30 GeV/c jet, compared to the reference. However, no significant broadening of 
the photon+jet azimuthal correlation is observed. 
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1 Introduction 

Parton scatterings with large momentum transfer produce energetic particles which can be 
used as "probes" to study the strongly interacting medium created in high-energy heavy ion 
collisions |TJ[2]|. The production of high transverse momentum (pj) partons and photons in 
"hard" processes occurs over very short time scales, r ~ 1/pj < O.lfm/c, and thus their 
yields can be potentially modified by final-state interactions occurring while they traverse the 
medium. Since the production cross sections of these energetic particles are calculable using 
perturbative quantum chromodynamics, they have long been recognised as particularly useful 
"tomographic" probes of the created medium (3HI1. 

In PbPb collisions at the Large Hadron Collider (LHC), the effects of the produced medium 
have been studied using back-to-back dijets which were observed to be significantly unbal- 
anced in their transverse momenta ffT0l - [T2| . The advantage of the large yield of dijets (as com- 
pared to photon+jet pairs) is, however, offset by a loss of information about the initial prop- 
erties of the probes, i.e. prior to their interactions with the medium. Correlating two probes 
that both undergo energy loss also induces a selection bias towards scatterings occurring at, 
and oriented tangential to, the surface of the medium. At leading order (LO), photons are 
produced back-to-back with an associated parton (jet) having close to the same transverse 
momentum. Furthermore, these photons do not strongly interact with the medium. The 
yields of isolated photons in PbPb collisions were found to match the expectation based on 
pp data and the number of nucleon-nucleon collisions, with a modification factor of Raa = 
0.99 ± 0.31(stat.) ± 0.26(syst.) [fT3l . Therefore, photon+jet production has been hailed as the 
"golden channel" to investigate energy loss of partons in the medium |T4l[T5ll . 

"Prompt photons" are photons produced directly in the hard sub-processes. Experimentally, 
events with enriched production of prompt photons are selected using an isolation require- 
ment, namely that the additional energy in a cone of fixed radius around the direction of the 
reconstructed photon be less than a specified value I1T31 . This restriction yields "isolated pho- 
tons" (7), which consist mostly of prompt photons produced directly in the initial hard scat- 
tering. Background photons from the decays of neutral mesons, such as n°, rj, and to, are 
suppressed by this isolation requirement, as they are predominantly produced via jet fragmen- 
tation. 

This Letter describes the first study of the jet energy loss using isolated-photon+jet pairs from 
PbPb data at a nucleon-nucleon centre-of-mass energy = 2.76 TeV. An integrated PbPb 

luminosity of f Cdt = 150 yb~ x was collected by the Compact Muon Solenoid (CMS) ex- 
periment during the 2011 running of the LHC. For comparison, a pp reference dataset with 
/ Cdt « 200 nb _1 at = 2.76 TeV was obtained in 2011. 

The goal of this analysis is to characterise possible modifications of jet properties as a function 
of centrality using isolated-photon+jet events in PbPb collisions. The properties of isolated- 
photon+jet pairs are studied via the azimuthal angular correlation in A(pj 7 = \(ft et — </> 7 | and 
the transverse momentum ratio given by Xj 7 = PjV pZ- Photons with transverse momentum 
of pj > 60GeV/c ar e sel ecte d in a pseudorapidity range of |// 7 | < 1.44, using isolation criteria 

These photons are then correlated with jets having p^ > 



detailed in Sections 



2.2 



and 



2.3 



30GeV/cand \rf et \ < 1.6. Parton energy loss due to induced gluon radiation can lead to a 
shift of the Xj 7 distribution towards lower values. In addition, parton energy loss can cause 
reconstructed jets to fall below the p j 1 > 30 GeV/ c threshold, leading to a reduction of the 
fraction of photons with an associated jet. 



Section 2 of this Letter begins with a description of the experimental setup as well as the event 
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triggering, selection and characterisation. The Monte Carlo simulation, the photon and jet 
reconstruction, and the analysis procedure are also described. The results and their systematic 
uncertainties are presented in Section 3, followed by a summary in Section 4. 

2 The CMS detector 

Particles produced in pp and PbPb collisions are studied using the CMS detector lTl6| . The 
central tracking system is comprised of silicon pixel and strip detectors that allow for the 
reconstruction of charged-particle trajectories in the pseudorapidity range \t]\ < 2.5, where 
rj = — ln[tan(0/2)] and 8 is the polar angle relative to the counterclockwise beam direction. 
Photons are reconstructed using the energy deposited in the barrel region of the PbWC>4 crystal 
electromagnetic calorimeter (ECAL), which covers a pseudorapidity range of < 1.479, and 
has a finely segmented granularity of Arj x A(p = 0.0174 x 0.0174. The brass/ scintillator hadron 
calorimeter (HCAL) barrel region covers \t]\ < 1.74, and has a segmentation of At] x Acp = 
0.087 x 0.087. Endcap regions of the HCAL and ECAL extend the |//| coverage out to about 3. 
The calorimeters and tracking systems are located within the 3.8 T magnetic field of the super- 
conducting solenoid. In addition to the barrel and endcap detectors, CMS includes hadron 
forward (HF) steel /quartz-fibre Cherenkov calorimeters, which cover the forward rapidity of 
2.9 < \r] | < 5.2 and are used to determine the degree of overlap ("centrality") of the two collid- 
ing Pb nuclei |I7|. A set of scintillator tiles, the beam scintillator counters, is mounted on the 
inner side of each HF for triggering and beam-halo rejection for both pp and PbPb collisions. 

2.1 Trigger and event selection 

Collision events containing high-px photon candidates are selected online by the CMS two- 
level trigger system consisting of the Level-1 (LI) and High Level Trigger (HLT). First, events 
are selected using an inclusive single-photon-candidate LI trigger with a transverse momen- 
tum threshold of 5GeV/c. Then, more refined photon candidates are reconstructed in the HLT 
using a clustering algorithm (identical to that used for offline analysis) applied to energy de- 
posits in the ECAL. Events containing a reconstructed photon candidate with p 2 > 40 GeV/ c are 
stored for further analysis. This HLT selection is fully efficient for events containing a photon 
with p2 > 50 GeV/ c and the analysis presented here includes all photons with p 2 > 60 GeV/ c. 

In order to select a pure sample of inelastic hadronic PbPb collisions for analysis, further offline 
selections were applied to the triggered event sample similar to ffTTl . Notably among these 
include requiring a reconstructed event vertex, and requiring at least 3 calorimeter towers in 
the HF on both sides of the interaction point with at least 3 GeV total deposited energy in 
each tower. Beam halo events were vetoed based on the timing of the +z and — z BSC signals. 
Additionally, events containing HCAL noise [ 18 1 are rejected to remove possible contamination 
of the jet sample. Details about this event selection scheme can be found in [10J. The number 
of events removed by these criteria are shown in Table [l] Analysis of the Monte Carlo (MC) 
reference, described in Section |2.2[ uses identical event selection, except for the calorimeter 
noise rejection, which is a purely experimental effect. 

The online trigger scheme for the pp data at 2.76 TeV is the same as that used for the CMS 
pp prompt photon analysis at 7 TeV |[T9|. The pp trigger requires at least one reconstructed 
electromagnetic cluster with a minimum transverse energy of 15 GeV/c. The offline criterion 
applied to select pp hadronic collision events is similar to previous CMS pp papers II20II . Apart 
from the trigger and hadronic collision selection the pp analysis uses the same event selections 
as the PbPb analysis [fl3l . 
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Table 1: The impact of the various event-selection criteria on the J Cdt = 150 mb~ PbPb data 
sample. In the third column, the percentages are with respect to the line above. The selections 
are applied in the sequence listed. Recall that A(pj 7 = \$ et — <p 7 \- 



Selection 


Events remaining 


% of previous 


Collision events with a photon of Pj > 40 GeV/c 


252576 




HCAL cleaning 


252317 


96.76 


Isolated photon candidate p\ > 60GeV/c, \tj\ < 1.44 


2974 


1.18 


Jet candidate p ] f > 30GeV/c, \t]\ < 1.6 


2198 


73.91 


A( Ph > \ n 


1535 


69.84 



For the analysis of PbPb events, it is important to determine the degree of overlap between 
the two colliding nuclei, termed collision centrality. Centrality is determined using the sum 
of transverse energy reconstructed in the HE The distribution of this total energy is used to 
divide the event sample into equal percentiles of the total nucleus-nucleus interaction cross 
section. These finer centrality bins are then combined into four groups; one containing the 10% 
most central events (i.e. those which have the smallest impact parameter of the two colliding 
Pb nuclei and which produce the highest HF energy); two encompassing the next most cen- 
tral 10-30% and 30-50% of the events; and finally one with the remaining 50-100% peripheral 
events. Centrality can also be characterised using the number of nucleons participating in the 
interaction, N par t (with N par t = 2 for pp). The corresponding N part values for a given centrality 
range are determined from a Glauber calculation [21 j. Detector effects are accounted for using 
a GEANT4 simulation [22] of events generated with a multi-phase transport model (AMPT) [23 1. 
A detailed description of the centrality determination procedure can be found in [fTOl . 

2.2 Monte Carlo simulation 

The production of high-px photons by LO processes and parton radiation and fragmentation 
channels with a high-px photon in the final state are simulated with PYTHIA fl24"| (version 6.422, 
tune Z2). Tune Z2 is identical to the Zl tune described in [25|, except that Z2 uses the CTEQ6L 
PDF while Zl uses CTEQ5L, and the cut-off for multiple parton interactions, p j_o, at the nomi- 
nal energy of = 1.8 TeV is decreased by 0.1 GeV/c. Modifications to account for the isospin 
effect of the colliding nuclei, i.e. the correct cross section weighting of pp, pn, and nn subcol- 
lisions |l26l| , is used. Events containing isolated photons are selected using the generator-level 
information of the PYTHIA events. The isolation criterion requires that the total energy within 
a cone of radius AR = y/ (A//) 2 + (A<^) 2 = 0.4 surrounding the photon direction be less than 
5 GeV. This selection is found to be equivalent to the experimental requirements for isolated 



photons described in Section 2.3 These events are then processed through the full CMS detec- 
tor simulation chain using the Geant4 package. In order to model the effect of the underlying 
PbPb events, the PYTHIA photon events are embedded into background events generated us- 
ing hydjet (v 1.8) |26|. This version of hydjet is tuned to reproduce event properties such as 
charged hadron multiplicity, pj spectra, and elliptic flow measured as a function of centrality 
in PbPb collisions. 

2.3 Photon reconstruction and identification 

Photon candidates are reconstructed from clusters of energy deposited in the ECAL, following 
the method detailed in Ref. [13]. The selected photon candidates are restricted to be in the 
barrel region of the ECAL by requiring a pseudorapidity limit of |?/ 7 | < 1.44 and are also 
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required to have a transverse momentum of p j > 60 GeV/ c. In addition, photon candidates 
are dropped if they overlap with any electron tracks. Electrons are identified by matching 
tracks with reconstructed ECAL clusters and putting a cut on the ratio of calorimeter energy 
over track momentum. The separation of the photon and electron is required to be within a 
search window of |f/ 7 — ^ Track | < 0.02 and \(p 7 — <^ Track | < 0.15. Anomalous signals caused 
by the interaction of heavily-ionising particles directly with the silicon avalanche photodiodes 
used for the ECAL barrel readout are removed, again using the prescription of Ref. [13]. The 
reconstructed photon energy is corrected to account for the material in front of the ECAL and 
for electromagnetic shower containment. An additional correction is applied to the clustered 
energy in order to remove the effects from the PbPb underlying event (UE). The size of the 
combined correction is obtained from the isolated photon PYTHIA + HYDJET sample and varies 
from 2-10%, depending on centrality and photon pZ. The effect of the corrections on the energy 
scale is validated by an analysis of the reconstructed Z boson mass observed in Z — > e e + 
decays in PbPb data as a function of centrality. 

Since the dominant source of neutral mesons is jet fragmentation with associated hadrons, 
a first rejection of neutral mesons mimicking a high-px photon in the ECAL is done using 
the ratio of hadronic to electromagnetic energy, H/E. The H/E ratio is calculated using the 
energy depositions in the HCAL and the ECAL inside a cone of AR = 0.15 around the photon 
candidate direction |[T9| . Photon candidates with H/E < 0.1 are selected for this analysis. A 
correction for the contribution from the remaining short-lived neutral mesons is applied later. 

To determine if a photon candidate is isolated, the detector activity in a cone of radius AR = 0.4 
with respect to the centroid of the cluster is used. The UE-subtracted photon isolation variable 
SumIso UE_sub , which is the sum of transverse energy measured in three sub-detectors (ECAL, 
HCAL, Tracker) minus the expected contribution from the UE to each sub-detector, as de- 
scribed in fl3|, is used to further reject photon candidates originating from jets. The mean 
of SumIso UE_sub for fragmentation and decay photons is ~ 20 GeV, while the distributions of 
SumIso UE_sub for isolated photons are Gaussians centred around and having widths vary- 
ing from 3.5 GeV for peripheral collisions to 8.5 GeV for central collisions. Candidates with 
SumIso UE " sub smaller than 1 GeV are selected for further study. A tightened isolation crite- 
rion for data (as compared to the 5 GeV applied for the MC) is used in order to minimise the 
impact of random PbPb UE fluctuations. A downward fluctuation in the UE contribution to 
SumIso UE_sub can inadvertently allow a non-isolated photon candidate to pass the isolation cut. 
From the PYTHIA + HYDJET sample, the fraction of photons surviving this tightened selection is 
estimated to be 70-85%, depending on centrality and photon pj, and is found not to be depen- 
dent on the angular or momentum correlation with the associated jet. The relative efficiency of 
SumIso UE ~ sub < 1 GeV compared to SumIso UE ~ sub < 5 GeV ranges from 82% (0-10% central- 
ity) to 90% (50-100% centrality). At the same time, the photon purity, central to peripheral, is 
74-83% for the SumIso UE ~ sub < 1 GeV cut, compared to 52-62% for the SumIso UE_sub < 5 GeV 
cut. 

Photon purities in each centrality interval are estimated using a two-component fit of the shape 
of the electromagnetic shower in the ECAL, a™, defined as a modified second moment of the 
electromagnetic energy cluster distribution around its mean rj position: 

_2 _ LiMVi-V) 2 

U 1 x (1) 

w t = max ( 0,4.7 + In 

where E, and r\\ are the energy and position of the z'-th ECAL crystal in a group of 5 x 5 crystals 
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centred on the one with the highest energy, E is the total energy of the crystals in the calcula- 
tion, and f\ is the average rj weighted by W{ in the same group [19J. The discrimination is based 
only on the pseudorapidity (i.e. longitudinal) distribution of the shower, which is aligned with 
the magnetic field direction. As a result, showers with a wider distribution in the transverse 
plane, which can originate from photons converted to e + e pairs in the detector material, are 
not eliminated. The shape of the distribution for the signal is obtained from photon+jet 
PYTHIA + HYDJET samples for each pj and centrality bin. The shape of the background dis- 
tribution is extracted from data using a background-enriched set of photon candidates with 
10 < SumIso UE_sub < 20 GeV. The estimated photon purity (one minus the nonphoton con- 
tamination) is 74-83% for photon candidates, which are required to have cr^ < 0.01. 

2.4 Jet reconstruction 

Jets are reconstructed by clustering particles measured with a particle-flow (PF) algorithm l|27| , 
using the anti-Zcj sequential recombination algorithm with a distance parameter of R = 0.3 [28J. 
The jets used in the analysis are required to have p ^ f > 30 GeV/ c and | rf et \ < 1 .6 to ensure high 
reconstruction efficiency. Jets within R < 0.3 around a photon are removed in order not to cor- 
relate the photon with itself. Details of the jet reconstruction procedure and its performance can 
be found in fl2~|. The small value of R, compared to a more typical R = 0.5-0.7 used to analyse 
pp events, helps to minimise sensitivity to the UE contribution, and especially its fluctuations. 
The energy from the UE is subtracted using the same method as employed in lH0llT2"| and origi- 
nally described in [29 1. The jet energy resolution can be quantified using the Gaussian standard 
deviation a of p j eco / p j en , where pj eco is the UE-subtracted, detector-level jet energy, and p j en 
is the generator-level jet energy without any contributions from a PbPb UE. The magnitude 
of this resolution is determined using PYTHIA + HYDJET simulation propagated through the 
detector using Geant4. Compared to direct embedding into PbPb events, this method avoids 
uncertainties associated with the detector versus MC geometry alignment, which is especially 
difficult to achieve accurately with finely segmented pixel trackers. The UE produced by HYD- 
JET with Geant4 has been checked against the data by observing the energy collected inside 
randomly oriented cones with the same radius as the distance parameter in the jet algorithm. 
Data and MC are found to be well matched. The dependence of the jet resolution, a, defined as 
the standard deviation of the reconstructed over the event generator p ^ t , can be parametrised 
using the expression 

a = C © . © (2) 



where © indicates a sum in quadrature, and the quantities C, S, and N are fitted parameters 
(Table [2j. The first two terms of the parametrisation are determined from PYTHIA simulation, 
and the third term, which represents background fluctuations (not corrected for the flow direc- 
tion), is determined from PYTHIA + HYDJET simulation. 

Because the effects of the UE for jets found in PbPb events are subtracted, corrections to the 
mean reconstructed jet energy are derived from pp data and PYTHlA-only simulation (i.e. with- 
out HYDJET) [30]. Studies of the performance of jet reconstruction in PYTHIA + HYDJET events 
show that no additional centrality-dependent energy correction is needed. 

The jet reconstruction efficiency is defined as the fraction of simulated PYTHIA jets which are 
correctly reconstructed when embedded into a HYDJET event. The efficiency is found to be 
greater than 90% for jets within the selected pj and t] range for all centralities. For the analysis 
of the pp sample, the same PbPb jet reconstruction algorithm is used. The performance of the 
jet reconstruction in peripheral PbPb events is found to approach that for the pp simulation. 
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Table 2: Parameters of the functional form for the jet energy resolution a (pj eco / Pj en ) given in 
Eq. 0, obtained from Geant4 simulation of pythia pp jets and from pythia jets embedded 
in hydjet events for various PbPb centralities (indicated by the % ranges in parentheses). The 
units of S are V GeV/c and the units of N are GeV/c. 



c 


S 


N(pp) 


N (50-100%) 


N (30-50%) 


N (10-30%) 


N (0-10%) 


0.0246 


1.213 


0.001 


0.001 


3.88 


5.10 


5.23 



2.5 Analysis procedure 

To construct photon+jet pairs, the highest p\ isolated photon candidate in each selected event 
is associated with every jet in the same event. The photon+jet pairs constructed in this way 
contain background contributions that need to be subtracted before using them to study energy 
loss effects on the jet produced in the same scattering as the photon. The dominant background 
contributions are photons from meson decays which pass the isolation requirement and the 
combinatoric background where the leading photon is paired with a jet not originating from 
the same hard scattering. The combinatoric background includes misidentified jets which arise 
from fluctuations of the underlying event as well as real jets from multiple hard interactions in 
the collision. 

The background contributions from decay photon and fake jets are estimated separately with 
methods that are data-driven and are subtracted from the photon+jet pair sample. 

The estimation of the yield and the kinematic characteristics of decay photons contained in 
the isolated-photon sample is based on the shower shape distributions for the analysed ECAL 
clusters. The ECAL clusters originating from high-px meson decays correspond to two photons 
that are reconstructed as a single wide cluster. Events with a large shower width (0.011 < cr„„ < 
0.017, see Eq. ^ are used to determine the contributions of the decay photon background to 
the Acpjy and xj 7 observables. The background shape obtained from this procedure is scaled 
according to the background-photon fraction, which is estimated from a fit of the shower shape 
distribution. The estimated background contribution fraction (which is equal to 1 — purity) is 
then subtracted from the yield for the signal events, which have a small shower width (cr„^ < 
0.01). 

The background contribution due to photon+jet pairs arising from fake jets or multiple hard 
scatterings is also subtracted. It is estimated by correlating each isolated highest-px photon 
from the triggered photon+jet sample to jets found in a different event selected randomly from 
a set of minimum bias PbPb data. The random event used in the pairing is chosen to have the 
same centrality as the photon+jet candidate event. The fake jet background estimated in this 
way has a flat distribution in A(pj 7 . The effect of this background is biggest in the most central 
events where, on average, approximately 20% of the jets paired with each photon candidate 
are estimated to be fake jets. The estimated distributions of A(/>j 7 and *j 7 for photons paired 
with fake jets, found using this random pairing of events, are subtracted from the distributions 
coming from the same-event photon+jet sample to obtain the final results. 

To determine the sensitivity to a potentially modified jet fragmentation, which may cause the 
reconstructed jet energy scale to deviate from the pythia derived calibration, the MC studies 
were repeated but now using PYQUEN 11261 jets embedded into hydjet. pyquen simulates par- 
ton energy loss by radiative and collisional mechanisms, where a portion of the original parton 
energy is redistributed into gluons that are found largely outside the cone of the surviving jet. 
The PYQUEN+HYDJET events were run through the full detector simulation and then recon- 
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structed with the standard analysis. The jet modification in PYQUEN produces a photon+jet 
momentum imbalance comparable to that observed in our measurement (although in detail, 
with different xj 7 distribution and N par t dependence). The extracted momentum imbalance 
was found to reproduce the generator level imbalance well within the statistical uncertainties. 
We also note that a similar insensitivity to differences among QCD fragmentation was found 
previously by studying the jet energy scale from separate pythia gluon and light quark jets, 
which differ significantly in their fragmentation patterns [31 1. The standard analysis using PF 
jets was cross checked using jets reconstructed with only information from the calorimeters. 
This alternative analysis has different corrections for jet energy scale and resolution and a dif- 
ferent sensitivity to low momentum tracks. The two analyses give comparable results for the 
photon+jet observables. 

3 Results 

3.1 Photon+jet azimuthal correlations 

Possible medium effects on the back-to-back alignment of the photon and recoiling jet can be 
studied using the distribution of the number of photon+jet pairs, Nj 7 , as a function of the 
relative azimuthal angle, A(pj 7 , normalised the total number of pairs, (Nj 7 )~ 1 dNj 7 /dA(pj 7 . 
Figure [l] shows distributions of A(pj 7 for PbPb data in four centrality bins, ranging from pe- 
ripheral events (50-100%, Fig.[lji) to the most central events (0-10%, Fig.j^l). The PbPb data 
are compared to PYTHIA + hydjet simulation and pp data. For both PbPb data and MC dis- 
tributions, the jet is found to be well aligned opposite to the photon direction, with a clear 
peak at A(pj 7 = u. The shape of the Acpj 7 correlation peak is similar in PbPb data and MC. 
The apparent excess in the tail of the 0-10% data was investigated and deemed statistically not 
significant compared to the subtracted background. To study the centrality evolution of the 
shape, the distributions are fitted to a normalised exponential function: 

1 dN h e (&<P-n)/<T 
N^ 7 dA(p }7 ~ (l-e- n/a )a 

The fit is restricted to the exponentially falling region Acp > 2n/3. The results of this fit for 
PbPb data are shown in Fig. |2j where the width of the azimuthal correlation (c in Eq. (J3}, de- 
noted cr(A(f)j 7 ) in Fig. [2j is plotted as a function of centrality and compared to pp and PYTHIA 
+ HYDJET fit results. The resulting cr(Acpj 7 ) values in PbPb do not show a significant centrality 
dependence within the present statistical and systematic uncertainties. For central PbPb colli- 
sions, cr(A(pj 7 ) is similar to the PYTHIA reference based on the Z2 tune, and comparison with 
other PYTHIA tunes shows a theoretical uncertainty that is larger than the difference between 
the data and MC. Comparing the PYTHIA tune Z2 with tune D6T ||32[[33l shows an 8% difference 
in cr(A(pj 7 ), which is expected because these two tunes differ in their parton shower ordering 
resulting in a different A(p correlation. The large statistical uncertainty in the o~(A(pj 7 ) extracted 
from the pp data at 2.76 TeV does not allow a discrimination between these two PYTHIA tunes. 
Both the Z2 and D6T tunes matched the shape of the azimuthal dijet correlation measured in 
pp collisions at 7 TeV [34] at about the 10% level in the region Acp > 2zr/3. The result that 
cr(A(pj 7 ) is not found to be significantly modified by the medium is consistent with the earlier 
observation of an unmodified A(p correlation in dijet events [10J. 

3.2 Photon+jet momentum imbalance 

The asymmetry ratio Xj 7 = PjV Pj is used to quantify the photon+jet momentum imbalance. 
In addition to the jet and photon selections used in the A(pj 7 study, we further impose a strict 
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Figure 1: Azimuthal correlation A<£j 7 between the photon and associated jet after background 
subtraction. The area of each distribution is normalised to unity. All panels show PbPb data 
(filled circles) compared to pp data at 2.76 TeV (filled squares), and to the pythia + hydjet MC 
simulation (shaded histogram) in bins of increasing centrality left to right. The error bars on 
the points represent the statistical uncertainty. 

A<pj 7 > g 7r cut to suppress contributions from background jets. Note that photon+jet pairs for 
which the associated jet falls below the 30 GeV/ c threshold are not included in the xj 7 calcu- 
lation. This limits the bulk of the xj 7 distribution to xj 7 > 0.5. Figure [3] shows the centrality 
dependence of xj 7 for PbPb collisions as well as that for pythia + hydjet simulation where 
PYTHIA contains inclusive isolated photon processes. The (xj 7 ) obtained from PYTHIA tunes Z2 
and D6T agree to better than 1%. Overlaid in the peripheral bin is the (xj 7 ) for 2.76 TeV pp data, 
showing consistency to the MC reference. However the poor statistics of the pp data and the 
50-100% PbPb centrality bin do not offer a strong constraint on a specific MC reference. How- 
ever, further studies using the 7 TeV high statistics pp data showed a good agreement in (xj 7 ) 
between data and PYTHIA, justifying the use of pythia + hydjet as an un-modified reference. 
The dominant source of systematic uncertainty in (x/ 7 ) is the relative photon+jet energy scale. 
Its impact on the probability density of xj 7 is approximately 10% for the intermediate region 
of 0.6 < Xjy < 1.2. The normalisation to unity causes a point-to-point anticorrelation in the 
systematic uncertainties, where the upward movement of the probability density at small xj 7 
has to be offset by the corresponding downward movement at large Xj 7 . This is represented by 
the separate open and shaded red systematic uncertainty boxes in Fig. [3] For a given change in 
the energy scale, all points would move together in the direction of either the open or shaded 
red box. The N par t dependence of the mean value (xj 7 ) is shown in Fig.|l|a). 

While the photon+jet momentum ratio in the PYTHIA + HYDJET simulation shows almost no 
change in the peak location and only a modest broadening, even in the most central PbPb 
events, the PbPb collision data exhibit a change in shape, shifting the distribution towards 
lower xj y as a function of centrality. It is important to note that, as discussed above, the limita- 
tion of Xj 7 > 0.5 limits the degree to which this distribution can shift. 

3.3 Jet energy loss 

To study the quantitative centrality evolution of the energy loss, the average ratio of the jet and 
photon transverse momenta, (xj 7 ), is shown in Fig. |4ja). While the photon+jet mean momen- 
tum ratio in the PYTHIA + HYDJET simulation exhibits a roughly centrality-independent value 
of (x h ) = 0.847 ± 0.004(stat.) - 0.859 ± 0.005(stat), the ratio is (x ]7 ) = 0.73 ± 0.02(stat.) ± 
0.04(syst.) in the most central PbPb data, indicating that the presence of the medium results in 
more unbalanced photon+jet pairs. 

It is important to keep in mind that the average energy loss of the selected photon+jet pairs does 
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Figure 2: Fitted Ac/>/ 7 width (a in Eq. (|3]l) between the photon and associated jet after back- 
ground subtraction as a function of N par t- The fit range was restricted to Ac/>/ 7 > %n. The 
yellow boxes indicate point-to-point systematic uncertainties and the error bars denote the sta- 
tistical uncertainty. 
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Figure 3: Ratio of p T between the photon (pj > 60 GeV/c) and jet (p'j > 30 GeV/c, A<pj 7 > 
g7r) after subtracting background. The area of each distribution is normalised to unity. All 
panels show PbPb data (filled circles) compared to pp data at 2.76 TeV (filled squares), and 
to the PYTHIA + HYDJET MC simulation (shaded histogram) in bins of increasing centrality 
left to right. The error bars on the points represent the statistical uncertainty. See text for an 
explanation of the open and shaded red systematic uncertainty boxes. 
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3 Results 



not constitute the full picture. There are genuine photon+jet events which do not contribute to 
the (xjy) distribution because the associated jet falls below the p^ > 30GeV/c threshold. To 
quantify this effect, Fig. Qb) shows Rj 7 , the fraction of isolated photons that have an associated 
jet passing the analysis selection. The value of Rj 7 is found to decrease, from Rj 7 = 0.685 ± 
0.008(stat.)-0.698 ± 0.006(stat.) for the pythia + hydjet reference, as well as pp and peripheral 
PbPb data, to the significantly lower R }7 = 0.49 ± 0.03(stat.) ± 0.02(syst.)-0.54 ± 0.05(stat.) ± 
0.02(syst.) for the three PbPb bins above 50% centrality. 

An analysis with a lower pj cutoff on the associated jet energy would result in values of Rj 7 
closer to unity. This would shift the cutoff at low xj 7 in Fig. 3 closer to zero. It is likely, although 
not certain, that these additional events would result in a larger deviation in Xj 7 between the 
PbPb data and the reference shown in Fig. 4(a). 



3.4 Systematic uncertainties 

Photon purity, reconstruction efficiency, and isolation, as well as the contamination from e and 
fake jets contribute to the systematic uncertainties of the photon+jet azimuthal correlation and 
the observables related to momentum asymmetry, {xj 7 ) and Rj 7 . Additionally, the momentum 
asymmetry observables are also influenced by the relative photon and jet energy calibrations. 
For the measurement of c(A</>), the uncertainty due to the photon angular resolution is negli- 
gible, less than 10~ 5 . 
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Figure 4: (a) Average ratio of jet transverse momentum to photon transverse momentum, 
(xj 7 ), as a function of N par t- The empty box at the far right indicates the correlated systematic 
uncertainty, (b) Average fraction of isolated photons with an associated jet above 30 GeV/ c, 
Rj 7 , as a function of N par t- In both panels, the yellow boxes indicate point-to-point systematic 
uncertainties and the error bars denote the statistical uncertainty. 



The uncertainty in the relative photon+jet energy scale consists of four main contributions. The 
first one comes from the 2% relative uncertainty of the jet energy scale in the barrel for 30 < 
Pj* < 200GeV/c, when compared with the ECAL energy scale II30II . The second contribution 
is the residual data-to-MC energy scale difference in pp collisions, which is not corrected for in 
this analysis, for which we quote the 2% maximum relative uncertainty which applies in the 
range \rf et \ < 1.6. Thirdly, the additional uncertainty for the jet energy scale in the presence 
of the UE is determined to be 3% for the 30 to 100% and 4% for the to 30% centrality range, 
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using the embedding of PYTHIA isolated photon+jet pairs into HYDJET. The fourth contribution 
is the effect of heavy ion background on the ECAL energy scale, which is determined from 
Z — > e~e + mass reconstruction, after applying the PbPb ECAL correction. This results in a 
relative uncertainty of 1.5%, comparable to the pp uncertainty (obtained via rc° and tj — > 77). 

The absolute photon energy scale uncertainty, estimated to be 1.5% using Z decays as described 
above, will also affect the threshold of our photon kinematic selection. Similarly, the lower 
transverse momentum cutoff for jets is sensitive to their absolute energy scale. For CMS, the 
energy of jets is calibrated by measuring the relative photon+jet energy scale in pp collisions, 
and therefore the uncertainty in jet energies is the quadrature sum of the uncertainties in the 
relative jet-to-photon energy scale and the absolute photon energy scale. 

The uncertainty of the photon purity measurement using the cr™ template fitting is estimated 
by (a) varying the selection of sideband regions that is used to obtain the background template 
and (b) shifting the template to measure the signal template uncertainty. These result in an 
estimated uncertainty on the photon purity of 12% and 2%, respectively. Systematic effects due 
to photon reconstruction efficiency are estimated by correcting the data using the efficiency de- 
rived from the MC simulation, and comparing the result with the uncorrected distribution. The 
contribution of non-isolated photons (mostly from jet fragmentation) that are incorrectly deter- 
mined to be isolated in the detector due to UE energy fluctuations or detector resolution effects 
is estimated using PYTHIA + hydjet simulation. The difference of photon+jet observables ob- 
tained from generator level isolated photons and detector level isolated photons is taken to be 
the systematic uncertainty resulting from the experimental criterion for an isolated photon. 

The current analysis removes contamination from fake jets purely by subtracting the back- 
ground estimated from event mixing. A cross-check of this subtraction has been performed 
using a direct rejection of fake jets via a fake jet discriminant. The discriminant sums the p\ 
of the jet core within R < 0.1 around the jet axis and determines the likelihood that the recon- 
structed jet is not the result of a background fluctuation. Both techniques for fake jet removal 
agree within 1% for the observables studied. The effect of inefficiencies in the jet finding is 
estimated by repeating the analysis and weighting each jet with the inverse of the jet finding 
efficiency as a function of \ 

Table 3: Relative systematic uncertainties for a(A(pj 7 ) for pp data and each of the PbPb central- 
ity bins. 



Source 


PP 


PbPb 50-100% 


PbPb 30-50% 


PbPb 10-30% 


PbPb 0-10% 


7 px threshold 


3.0% 


3.0% 


3.0% 


2.0% 


1.2% 


Jet pj threshold 


1.3% 


1.3% 


0.2% 


0.5% 


2.4% 


7 efficiency 


0.8% 


0.8% 


0.3% 


0.3% 


0.3% 


Jet efficiency 


0.6% 


0.6% 


0.7% 


0.4% 


0.3% 


Isolated 7 definition 


0.7% 


0.7% 


1.6% 


2.0% 


0.5% 


7 purity 


6.8% 


6.8% 


2.7% 


0.5% 


0.9% 


e~,e + contamination 


0.5% 


0.5% 


0.5% 


0.5% 


0.5% 


Fake jet contamination 


0.3% 


0.3% 


0.1% 


0.2% 


1.2% 


Jet <p resolution 


0.5% 


0.5% 


0.5% 


0.5% 


0.5% 


cr fitting 


0.3% 


0.3% 


0.1% 


0.1% 


0.1% 


Total 


7.7% 


7.7% 


4.5% 


3.0% 


3.2% 



Tables Big andg summarise the relative systematic uncertainties for a(A<p), {xjy), and Rj 7 , 
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3 Results 



Table 4: Relative systematic uncertainties for (x/ 7 ) for pp data and each of the PbPb centrality 
bins. The uncertainties due to the pp 7-jet relative energy scale and 7 purity are common to all 
of the measurements and are quoted as a correlated uncertainty. 



Source 


PP 


PbPb 50-100% 


PbPb 30-50% 


PbPb 10-30% 


PbPb 0-10% 


7-jet rel. energy scale 


2.8% 


4.1% 


5.4% 


5.0% 


4.9% 


7 px threshold 


0.6% 


0.6% 


0.6% 


0.6% 


1.3% 


Jet px threshold 


0.7% 


0.7% 


1.9% 


1.9% 


2.0% 


7 efficiency 


< 0.1% 


< 0.1% 


< 0.1% 


0.1% 


0.2% 


Jet efficiency 


0.5% 


0.5% 


0.6% 


0.6% 


0.5% 


Isolated 7 definition 


0.1% 


0.1% 


0.7% 


0.4% 


2.0% 


7 purity 


2.2% 


2.2% 


1.9% 


2.4% 


2.7% 


e~,e + contamination 


0.5% 


0.5% 


0.5% 


0.5% 


0.5% 


Fake jet contamination 


0.1% 


0.1% 


0.1% 


0.2% 


0.1% 


Total 


3.7% 


4.8% 


6.2% 


6.0% 


6.4% 


Correlated (abs., rel.) 


3.6% 


3.6% 


3.6% 


3.6% 


3.6% 


Point-to-point 


0.9% 


3.2% 


5.1% 


4.8% 


5.3% 



Table 5: Relative systematic uncertainties for the fraction of photons matched with jets, Rj 7 , for 
pp data and each of the PbPb centrality bins. 



Source 


PP 


PbPb 50-100% 


PbPb 30-50% 


PbPb 10-30% 


PbPb 0-10% 


7 px threshold 


2.0% 


2.0% 


1.9% 


1.3% 


2.1% 


Jet px threshold 


1.4% 


1.4% 


2.3% 


2.6% 


2.7% 


7 efficiency 


0.2% 


0.2% 


0.2% 


0.5% 


0.5% 


Jet efficiency 


1.5% 


1.5% 


1.7% 


1.8% 


2.1% 


Isolated 7 definition 


0.2% 


0.2% 


0.6% 


1.3% 


0.8% 


7 purity 


2.3% 


2.3% 


1.9% 


0.2% 


0.9% 


e~,e + contamination 


0.5% 


0.5% 


0.5% 


0.5% 


0.5% 


Fake jet contamination 


0.4% 


0.4% 


0.8% 


1.0% 


1.4% 


Total 


3.7% 


3.7% 


4.1% 


3.9% 


4.5% 
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respectively, for the pp data and for each of the PbPb centrality bins used in the analysis. For 
(xj 7 ), the uncertainties are separated into a correlated component that is common to all cen- 
trality bins and a component that represents the point-to-point systematic uncertainty. The 
common correlated uncertainty is obtained by combining the pp jet energy scale uncertainty 
with the photon purity uncertainty. This absolute uncertainty of 3.6% was used as the corre- 
lated uncertainty for all PbPb centrality bins. 

4 Conclusions 

The first study of isolated-photon+jet correlations in PbPb collisions at ^/snn = 2.76 TeV has 
been performed as a function of collision centrality using a dataset corresponding to an inte- 
grated luminosity of 150 pb -1 . Isolated photons with pj > 60GeV/c were correlated with jets 
with > 30GeV/c to determine the width of the angular correlation function, c(A<pj 7 ), the 
jet/photon transverse momentum ratio, Xj 7 = p^V p\, and the fraction of photons with an 
associated jet, Rj 7 . The PbPb data were compared to both pp data and a pythia + hydjet 
MC reference which included the effect of the underlying PbPb event but no parton energy 
loss. No angular broadening was observed beyond that seen in the pp data and MC refer- 
ence at all centralities. The average transverse momentum ratio for the most central events 
was found to be (*/ 7 )o-io% = 0-73 ± 0.02(stat.) ± 0.04(syst.). This is lower than the value of 
0.86 seen in the pp data and predicted by pythia + hydjet at the same centrality. In addi- 
tion to the shift in momentum balance, it was found that, in central PbPb data, only a fraction 
equal to R/ 7 = 0.49 ± 0.03 (stat.) ± 0.02 (syst.) of photons are matched with an associated jet at 
A(pj 7 > \tc, compared to a value of 0.69 seen in pythia + hydjet simulation. Due to the hot 
and dense medium created in central PbPb collisions, the energy loss of the associated parton 
causes the corresponding reconstructed jet to fall below the p^ 1 > 30GeV/c threshold for an 
additional 20% of the selected photons. 
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